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Abstract. We present a maximum likelihood analysis of cosmological parameters from measurements of the 
aperture mass up to 35 arcmin, using simulated and real cosmic shear data. A four-dimensional parameter space 
is explored which examines the mean density £Im, the mass power spectrum normalization as, the shape parameter 
r and the redshift of the sources z s . Constraints on Qm and as (resp. F and z s ) are then given by marginalizing 
over r and z 3 (resp. Qm and as)- For a flat ACDM cosmologies, using a photometric redshift prior for the 
sources and F G [0.1,0.4], we find a s = (0.57 ± 0.04) n^ 24=F ° 18)nM " ' 49 at the 68% confidence level (the error 
budget includes statistical noise, full cosmic variance and residual systematic). The estimate of F, marginalized 
over Qm £ [0.1,0.4], as € [0.7,1.3] and z 3 constrained by photometric redshifts, gives F = 0.25 ±0.13 at 68% 
confidence. Adopting h = 0.7, a flat universe, F = 0.2 and Q m = 0.3 we find crg=0.98 ±0.06 . Combined with 
CMB, our results suggest a non-zero cosmological constant and provide tight constraints on Qm and erg. We 
finaly compare our results to the cluster abundance ones, and discuss the possible discrepancy with the latest 
determinations of the cluster method. In particular we point out the actual limitations of the mass power spectrum 
prediction in the non-linear regime, and the importance for its improvement. 
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1. Introduction 

In the standard cosmological picture, the structures in the 
Universe grow from the gravitational collapse of initial 
Gaussian density perturbations. The properties of mass 
distribution at low redshift are expected to express the 
latest and one of the most explicit footprint of the forma- 
tion process, so their description from cosmological sur- 
veys can be challenged against theoretical predictions re- 
sulting from this paradigm. For example, a direct obser- 
vation of the mass distribution in structures is believed to 
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be an unequivocal test of cosmological scenario of struc- 
ture formation. If so, the weak gravitational lensing pro- 
duced on distant galaxies by large scale structures is a 
direct probe of dark matter, regardless the light distribu- 
tion. It is therefore a robust technique to challenge the 
current cosmological models. In particular, it can reliably 
probe small angular scale and look into the transition to 
the quasi-linear and non-linear regimes, where comparison 
between observations and cosmological models are still dif- 
ficult. 

The cosmological origin of the coherent distortion 
fields detected in cosmic shear surveys is now firmly es- 
tablished ([Bacon et al., 2000 ; Hacmmcrlc ct al ., 2001 



Hockstra ct al., 2002) ; |Kaiser et al., 2000| |Maoli et al 



de l'TTrivprs (TNSTT) nf the Hpnt.rp National Hp la Recherche 



2001: Pen ct al., 2001: Rhodes et al., 2001 



Scientii .que (CNRS) and the University of Hawaii fUH). and et al., 2000| , |2001b| ; |Wittman et al., 2000[) . |Van Waerbeke 
at the European Southern Observatory telescopes Very Large et al. (2001b| ) have shown that the measurements provided 



Van Waerbeke 



Telescope (VLT) and the New Technology Telescope (NTT) . by different statistical estimators of distortion signal are 
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consistent with the gravitational lensing hypothesis with a 
high confidence level, so that present-day data can already 
constrain cosmological parameters. Their joint estimate of 
the mass density J7m and the power spectrum normaliza- 
tion as led to consistent results with the cluster abundance 



constraints (Picrpaoli ct al., 2001) and confirmed earlier 



tentative obtained by Maoli et al. (2001) and Rhodes et al. 
(200l[j]using ESO-VLT/CFHT and HST data respectively. 
A recent measurement done on a shallow survey (there- 
fore very different in depth) confirmed also this agreement 
([Hockstra et al., 2001|, |2002|). 



So far, the cosmological parameter estimation from 
cosmic shear relied on prior knowledge of the slope of the 
mass power spectrum T and/or the mean redshift z s of the 
lensed galaxy population. In fact, the statistical proper- 
ties of cosmic shear depends significantly on these quanti- 



ties (Kaiser, 1992 



Bernardeau et al., 1997 ; Jain fc Seljak 



1997|), |so any prior on these parameters may have a se- 
rious impact on the cosmological parameter estimation. 
For instance changing the shape of the power spectrum 
in either direction would favor low or high matter densi- 
ties, by changing the normalization accordingly. This am- 
biguity expresses a degeneracy between the normalization 
and the mass density, which depends on the choice of T 
(|Van Waerbeke et al., 2001b|). [Jain fc Seljak (1997|) ad- 



dressed this issue by pointing out that a measurement of 
the cosmic shear in both linear and non-linear scales could 
break the degeneracy, so that one in principle recover si- 
multaneously r, us and CIm from the shear variance alone. 
Unfortunately, the redshift of the sources is also a strongly 
degenerate parameter with ag, which definitely hampers 
shear variance analysis to provide unequivocal discrimina- 
tion of cosmological models. In fact, stringent constraints 
on the cosmological parameters from the shear variance 
are possible only with an accurate knowledge of the source 
redshifts and a measurement which extends over a large 
range of scales. 

In this paper we carry out a full maximum likelihood 
analysis of cosmic shear data over the four parameters 
O^f, as, r, and z s for flat and open cosmologies. Using 
both simulations and observations, we study slices and 
projections in this parameter space and discuss the re- 
liability of cosmological constraints which are based on 
catalogues having similar size and depth as current cos- 
mic shear surveys. In particular we give an estimate of 
and (j 8 by marginalizing over the power spectrum shape 
and sources redshifts. The improvement of our knowledge 
of the source redshift is crucial to gain better accuracy 
for the other parameters. In this work we use photometric 
redshifts ^] to put priors on the source redshift distribution. 

The paper is organized as follows. Section 2 presents a 
brief summary of some theoretical concepts and introduces 
the shear estimators used throughout the paper. Section 3 
describes the data and how shear quantities were obtained 
from the survey catalogue. The likelihood method and the 
details of the priors are presented in Section 4. Section 5 



shows and discusses the results on the parameter estimates 
on both simulated and real surveys. Finally, conclusions 
are presented in Section 6. 

2. Theory 



Following the notation in Schneider et al. (1998| ) , we define 
the power spectrum of the convergence as 



P K {k) = 



4 "o 



a 2 (w) 



P: 



:\D 



Jk(w) 



; w 



, / / /n,/k(w -w) 
dw n(w ) — , jt 

JKW) 



(1) 



where /k(w) is the comoving angular diameter distance 
out to a distance w (wh is the horizon distance), and 
n(w(z)) is the redshift distribution of the sources given 
in Eq.(|12j). P3_o(fc) is t he non-linear mass power spectrum 
computed according to Peacock fc Dodds (1996 ) , and k is 
the 2-dimensional wave vector perpendicular to the line- 
of-sight. The top- hat shear variance (smoothing window 
of radius 9 C ) and the shear correlation function can be 
written as 



<7 2 > 



2 

(77), = ± 



dk kP K (k)J (kd). 



(2) 



(3) 



Because the weak distortion field can be generated by non- 
lensing mechanisms, it is important to measure separately 
the E and B components of the shear. These modes were 
introduced initiall y to test for the gravita tional origin of 
the lensing signal ( Crittenden et al., 2000 ) since a poten- 
tial gravitational field is expected to produce only curl-free 
shear patterns (E mode). Any measurable B mode should 
be interpreted as a measurement of residual systematic 
in the data (Point Spread Function correction, intrinsic 
alignment or anything else) and must be removed from 
the weak lensing signal prior to cosmological interpreta- 
tion of cosmic shear data. 

The extraction of both modes is not trivial. The E and 
B-modcs decompositions of the top-hat shear variance, 
and of the shear correlation function given in Eq.(|2p|) are 



only defined up to a integration constant (see Crittenden 
|ct al., 200C ; Pen et al., 2001). This constant depends on 
the extrapolated cosmic shear signal either at small (< 30 
arc-seconds) or large (> 1 degree) scales. These boundary 
conditions turn out to be a severe limitation which ham- 
pers reliable derivations of both modes from our present- 
day data because we do not cover yet very large angular 
scales and we still suffer from systematics on very small 
scales which are not well understood. As pointed out by 



Pen ct al. (2001), the only unambiguous E and B mode 
decomposition is carried out by the aperture mass, M ap : 



derived from other data sets 



M a 



d 2 K (0) U(0), 



(4) 
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where k(9) is the convergence field, and U(9) is the zero 
mass aperture window ( [Schneider et al., 199S ) : 



U(6) 



9 



91 



(5) 



This estimator was introduced in Kaiser et al. (1994) to 
study clusters of galaxies, but it also has an important 
potential for cosmic shear analysis ( Schneider et al., 1998j ) . 



(Mf p ) can be calculated directly from the shear 7 with- 
out the need for a mass reconstruction. It can be rewritten 
as a function of the shear if we express 7 = (7*, 7r) m the 
local frame of the line connecting the aperture center to 
the gala xy. M ap can therefore be expressed as function of 



7t only (Miralda-Escude, 1991; Kaiser, 1992) 



M, 



ap 



d 2 6 lt {6) Q{0), 



(6) 
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where the filter Q(9) is given from U(9): 



Q{0) 



d6' 9' U(9') - U{9) 



(7) 



Fig. 1. Above is the aperture mass statistic (M 2 p ) (E- 
mode) and below the aperture mass (M 2 p )± computed 
with galaxies rotated by 45 degrees (i3-mode). Error bars 
are 1 — a statistical errors. 



The aperture mass variance is related to the convergence 
power spectrum (Eqj^)by: 



288 
tt0* 



dk 



P K {k)[J A {k9 c )f 



(8) 



The -B-mode is obtained by replacing 74 with 7 r in Eq. (g) . 
Although this estimator is robust and does not depend on 
an unknown integration constant, it is less sensitive to 
cosmological parameters th an the top-hat variance or th e 
shear correlation functions ( Van Waerbeke et al., 2001b|) . 



3. Measurements 

We use the observations done within the VIRMOS- 
DESCART project I by the VIRMOS f\ imaging and 
spectroscopic survey. The data cover an effective area 
of 8.5 sq.deg. in the I-band, with a limiting magnitude 
m] AB = 24.5. Technical details of the data set are given 
in Van Waerbeke et al. (2001b). We applied a bright mag- 
nitude cut at mi — 21 in order to exclude the foreground 
objects from the source galaxies. The shape of the galaxies 
are measured and analyzed as described in this paper, to 



^WiW-jetfii) ■ e t (0j) 
id 

WiW 3 e r {0i) ■ e r (8j) 
hi 



(9) 



where r = \Oi—0j\, and (et, e r ) are the tangential and ra- 
dial ellipticities defined in the frame of the line connecting 
a pair of galaxies. 

From Eq.(|)|), we define £+(r) and £_(r) which are re- 
spectively the sum and the difference of the two correlation 
functions: 



C+M = &(r) +Crr(r); £_(r) = &(r) - &. r (r). 



(10) 



which tvc refer tor technical details. 



Both (r) and £_ (r) are computed from a summation of 
the correlation function defined in Eq.(10), while the E 
and B modes aperture mass are derived by integration of 
the correlation functions with an appropriate window (see 
|Crittcndcn et al. (2000| ) for general derivations and Pen 



The location of the i-th galaxy is given by 0i, the ellip- 
ticity by e(0j) = (ei, £2), and its weight W{. The ellipticity 
is an unbiased estimate of the shear j(9i). The quantity 
measured on the data are the binned tangential and radial 
shear correlation functions. They are given by a sum over 
galaxy ipaiis 9j) 



et al. (2001 ) for a practical application to our filter). 



The E mode aperture mass is 

p20 c 

(ML)=ir rdrW{r)i+{r)+Ti 



20,. 



rdrW(r)£-(r), (11) 



where W(r) and W(r) are given in Crittenden et al 
(2000D ^|. The B-mode is obtained by changing the sign 



2 http://terapix.iap.fr/DESCART 

3 http://www.astrsp-mrs.fr 



Useful exp ressions using similar formalism as this work 
can be found in [Schneider et al. (200l| ) 
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of the second term in Eq. ( |ll| ) (which is equivalent to the 
45 degrees rotation test, or else jt — ► 7r)- 

Figure |l| shows the E mode (top) and B mode (bot- 
tom) measured in our galaxy sample. Using the B-mode 
measurement, we found the source of the residual systcm- 
atics at 3 — 4' report ed by Van Waerbeke et al. (2001b ) 
and Pen et al. (2001 ): it was caused by the third order 



polynomial fit to the PSF, which produced wings at the 
edge of the CCDs. A second order fitting removed most of 
the unwanted B mode contribution without spoiling the E 
mode signal. As shown in Figure ^, the residual system- 
atics are consistent with zero up to 10 arc-minutes and 
remains flat over the whole angular scale explored by the 
data. Clearly, the signal is dominated by the E mode con- 
tribution at least up to 25 arc-minutes. This demonstrates 
that signal produced by intrinsic alignment of galaxies is 
not detected at this level. 



4. Parameter Estimation 



90% of the Iab < 24.5 galaxies. Hence, most galaxies used 
for determining the photo- z distribution of galaxies up to 
Iab < 24.5 have reliable J and K photometric measure- 
ments to secure a collection of redshift on a very large 
sample of galaxies, which covers a broad magnitude range 
homogeneously spread over the whole field. The presence 
of the lensing cluster (z = 0.306) in the field only affects 
the redshift range 0.26 < z < 0.36. These data have been 
removed from the sample and the redshift distribution in- 
terpolated in this redshift range. The magnification bias 
may also change the redshift distribution of galaxies inside 
the very center of the cl uster where the gravi tational de- 
pletion is significant (see Athrcya et al., 2002| ). We there- 
fore also removed the central part (R < 40 arc-second) of 
the cluster from the sample. Since this region is also the 
most contaminated by the brightest cluster members, the 
depletion itself turns out to have no impact of the galaxy 
selection criterion. 

The photometric rcdshifts (z p hot) were measured us- 



4.1. R edshift distribution of galaxies in the 



ing the fitting algorithm hyperz developed by Bolzonella 
et al. (200C). Each z p i lot is inferred by comparing the 



VIRMOS-DESCART data 

We estimated the redshift distribution of our catalogue 
from a com bination of the Hubble Deep Fields North and 
South data ( |Fernandez-Soto et al., 1999| ; |Chen et al., 199§|) 
and VLT observations of the cluster MS1008-1224. Both 
HDF and MS1008-1224 observations are much deeper 
than the cosmic shear sample of galaxies considered in 
this work, so that magnitude measurements and photo- 
metric redshift estimates up to Iab = 24.5 are based on 
reliable data with high signal-to-noise ratio. 

The VLT MS1008-1224 galaxy sample comprises 
deep UBVRI observations, carried out by the Science 
Verification Team (SVT) at ESO/VLT with the FO RS1 
and FORS2 instruments ( |Appcnzcllcr et al., 1998] ) and 
deep J and K data obtained at the ESO/NTT with SOFI 
(Program 66.A-0316(A); PI Mellier). The extension of 
early deep SVT observations to U band with FORS2 and 
more recently to near infrared with SOFI, which has simi- 
lar field of view as FORS (5.5 arc-minutes against 6.8 arc- 
minutes), allow us to considerably improve the accuracy 
of photometric-redshifts of foreground, cluster and back- 
ground galaxies over the whole field. As compared to HDF 
the VLT/NTT observations are not as deep, but they pro- 
vide a much larger sample of galaxies because they cover 
of field of view 15 times larger than HDF. In total, 920 
galaxies having Iab < 24.5 and UBVRIJK data have 
been added to HDF data. 

The deep UBVRI data a re described at the ESO siteQ 
and in Athrcya et al. (2002). A complete description of 



spectral energy distribution of galaxies, as sampled by 
their UBVRIJK photometric flux, to a set of spectral 
templates representative of common late and early type 
galaxies which are followed with look-back time accord- 
ing to Bruzual fc Chariot' s evolution models (GISSEL98; 
Bruzual fe Chariot, 1993). The validation of h yperz is 
discussed at length in Bolzonella et al. (2000|) . It has 



been conclusively gauged against spectroscopic redshift 
on MS1008-1224 data by |Athreya et al. (2002| ). Details 
on photometric redshift techniques can be found in those 
papers. 

The compiled photometric distribution is shown on 
Fig.^. For the purpose of marginalization we parameterize 
this distribution with the following normalized function: 



n(z) = 



-, exp 



(12) 



the new J and K band data will be presented elsewhere 
( pavazzi fc et al., 2002 ). In brief, the exposure times of 
NTT/SOFI J and K bands were 5h30 and 6h respectively. 
The completeness limits are J = 23. and K — 22. and 
both J and K complete samples encompasses more than 



where a = 2 and (3 = 1.2. For these values of a and /?, 
the mean redshift is z s ss 2.1 z s and the median redshift is 
w 1.9 z s . We allowed z s to vary from 0.24 to 0.64, which 
corresponds to a mean redshift varying from 0.5 to 1.32. 
These two extreme models are shown on Fig.|^: they are 
clearly conservative bounds on the redshift distribution in 
the data. The curve on FigJ| shows the best fit model, 
with z s = 0.44 (z s = 0.9). 

4.2. Maximum Likelihood 

The dominant cosmological parameters for the 2-point 
cosmic shear statistics are the mean mass density Qm, 
the power spectrum normalization og, the shape parame- 



ter r and the redshift of the sources (see Bernardeau et al 



19971 ; |Jain fc Seljak, 1997| ; |van Waerbeke et al., 1999|) Our 



http:/ /www.hq.eso.org/s 



e/utlsv 



parameter space has therefore four dimensions, but we 
truncate the exploration volume inside a realistic range 
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Fig. 2. The histogram shows the photometric redshift 
distribution from the joint VLT and HDF fields. The 
thick solid line is the theoretical n(z) from Eq.([i"2|) with 
z s = 0.44. The low and high redshift thin solid lines cor- 
respond to z s = 0.24 and z s = 0.64, our extreme redshift 
distributions used in this paper. 
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Fig. 3. Cross-correlation coefficient between three mea- 
surement scales for the statistical noise (above) and 
the cosmic variance (below). The three scales are 1.3' 
(squares), 10.3' (triangles) and 25.1' (diamonds). The cos- 
mic variance cross-correlation is obtained from ray-tracing 
simulations ( Jain et al. (2000| ) ) . 



defined as Q M S [0.1,1], a s € [0.3,1.6] and T G [0.05,0.7] 
with a sampling of 10 x 14 x 14. For the analysis of the 
VIRMOS-DESCART data, the redshift of sources is pa- 
rameterized by Eq. (|l2l) with z s £ [0.24,0.64] and a sam- 
pling of 9. For the simulations, the sources are placed at 
redshift unity, therefore in the maximum likelihood anal- 
ysis we assumed we knew the shape (Dirac distribution), 
but we allowed the redshift z s to vary between 0.7 and 
1.2 (sampling of 6). In fact we found that the real shape 
of the source distribution does not matter, but the agree- 
ment with the mean redshift does. This parameter range 
box (CImi I\ &8, z s ) defines what we call the default prior 
box. The model predictions are then interpolated with an 
oversampling seven times better in each dimension. 

Let di be the data vector [i.e. the aperture mass (Mf p ) 
for different scales 0,), and rrii(QMiO's,T,z s ) the model 
predictions. The likelihood function of the data is 



L = 



(2^e 



112. 



exp 



(di - rn, ; )C L (di - m"] 7 



the bottom panel in Figure [j]) to the error of the signal. 
Given that there is no clearly identified scheme to deal 
with the residual systematics yet, this appears to be the 
safest and most conservative attitude. The diagonal part 
of the bias correlation matrix C& is therefore given by 
the B mode signal, and the off-diagonal terms follow the 
same correlation properties as the E mode (the E and B 
covariance matrices for the statistical noise are actually 
identical) . 

The cosmic variance covariance matrix C s is trickier to 
estimate. Assuming the field is Gaussian is too simplistic, 
since the observed scales are within the non-linear and 
weakly non-linear regimes, so in principle a complete de- 
scription of non-Gaussian contributions to the error bud- 
get cannot be carried out without detailed cosmological 
simulations. In order to avoid this heavy procedure, we fo- 
cused instead on a simpler alternative based on non-linear 
perturbation theory. It was pointed out in Scoccimarro 



(13) et ah (1999) that, for the convergence power spectrum, 



where n = 10 is the number of scales and C is the 10 x 10 
covariance matrix, 



Cij = ((di - m i ) T (d j - rrij)}. 



(14) 



C can be decomposed as C = C„ + C s + C&, where C„ 
is the statistical noise, C s the cosmic variance covariance 
m atrix and Cb the residual bias. C n has been measured 
in Pen et al. (2001 ), so we just reproduce here its gen- 
eral b ehavior: the top panel on figure shows the cross- 



the ratio of the Gaussian to Non-Gaussian errors is almost 
independent of scale, and close to 1 for any cosmology. We 
investigated whether this statement could be also valid in 
real space, using three ray-tracing simulations for three 
different cosmological models ( Jain et al., 200C| ). Figure 
I shows this ratio for ACDM, rCDM and OCDM. For 
scales larger than 3 arcmin., it is indeed nearly indepen- 
dent of cosmology. At smaller scales the ACDM model has 
a larger cosmic variance, but this is not important, since 
below a few arcminutes statistical noise do minates (see 
Figure p|). T herefore, although the result in Scoccimarro 



correla tion coefficient for 2, 10, and 35 arc-min with the 
other scales. In order to account for residual systematics, 
we decided to add quadratically the residual B mode (see 



et al. (1999| ) is clearly not valid for very small scales, it 



is still weakly sensitive to cosmology. We then approxi- 
mated C s in the following way: we compute the Gaussian 
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Fig. 4. Ratio between the aperture mass dispersion ob - 
tained from ray-tracing simulations ( Jain et al. (2000| )) 
and the dispersion obtained from the Gaussian field hy- 
pothesis for a survey of similar size. The three curves 
correspond to a rCDM (squares), OCDM (triangles) and 
ACDM (diamonds) , showing that the ratio is little depen- 
dent of the cosmological model above 3'. paper. 



cosmic variance for each model, then we convert it to a 

nnn-nanssian rnsmir varianrp using the awragp Vnrtnsis 



111 Fig tlre|| The prnss-rnrrelatinn rneffiripnt. is taken frnm 
the ray-tracing simulations. The different scales are rather 



correlated, as shown on Figure [j (bottom panel). As we 



shall see in Figure |9|, even a wrong estimate of the cosmic 
variance by a factor of two has no consequences on our 
parameter estimate, given that the errors are dominated 
by C„ and C b . 



5. Applications 

We now apply the likelihood analysis to simulated sky 
images and to the VIRMOS-DESCART data. 



5.1. Mock catalogues 

The mock catalogues are generated from simulated sky 



images following the procedure described in ( Erben et al 
2001) a simulated catalogue of galaxies is first lenscd and 



then u=cd to generate a CCD image of the sky. But in- 
stead of having a constant shear amplitude on each field, 
the distortio n of galaxies is i ntroduced using ray-tracing 
simulations ( Jain et al., 200C ) 

As for real sky surveys, the mock catalogues contain 



the fol owing features 



— galaxy intrinsic shape fluctuations; 

— masks; 

— noise from galaxy shape measurements and systemat- 
ics from PSF corrections , 




Fig. 5. Histogram of the Point Spread Function 
anisotropy of the stars in the simulated images (solid line) 
and in the VIRMOS-DESCART survey (dashed line). 



and the simulated images reproduce similar observational 
conditions as real data (PSF anisotropy, limiting magni- 
tude, luminosity functions, galaxy and star number den- 
sities, intrinsic ellipticity...). The simulated galaxies are 
then analyzed exactly in the same way as real data, fo l- 



lowin g the procedure described in ( |Van Waerbeke et al. 
20001 , 



2001b) 



We used two ray-tracing simulations from Jain et al 
(20001 ): one is OCDM, as described in Section 3, and the 
other is a rCDM with T = 0.21 and tt = 1. For each 
simulation we produced 11 square degrees of simulated 
sky images containing roughly 30 galaxies per arcmin 2 , 
with a pixel size of 0.2 arcsec. Figure || compare the star 
anisotropy between the simulated fields (solid line) and the 
data (dashed line). The likelihood function is computed 
for 11760 models (10 x 14 x 14 x 6). Figure || shows the 
results for the maximum likelihood analysis of these two 
simulation sets. We clearly converge to the right cosmo- 
logical model, which validates our likelihood approach for 
the data (see Section 5.2). However, we should point out 
here that the likelihood method assumes that our theoret- 
ical predictions are accurate compared to the precision of 
the measurements. Our simulations shows this is unfortu- 
nately not necessarily the case with today's lensing data 
sets. For instance, it was shown in Van Waerbeke et al 



(2001a) (Figure 2), that the non-linear predictions fails 
badly for the aperture mass with a rCDM model. This 
failure should not be a surprise: it was already noticed in 
the projected power spectrum in Jain et al. (2000|) (Figure 



and even the VIRGO simulations (see Jenkins et al 



1998, Figure 7) already noticed a mismatch between the 



3D non-linear predictions and the measured power spec- 
trum. In the case of our rCDM simulation, the potential 
problem is an overestimate of the power spectrum nor- 
malization erg. This is illustrated in Figure fj] where we 
compare the measured power to the Peacock & Dodds 
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Fig. 6. Constraints on as and Q,m for the OCDM (left) and tCDM (right) simulations. The gray levels and the 
dashed contours show the -simulation- default prior constraints, with confidence levels of 0.68, 0.95 and 0.999. The 
true model, indicated by a dark triangle, is ~ 0.3, as — 0.85 (left) and f2 = 1, erg = 0.6 (right). The thick solid line 
contours are for a prior T £ [0.1, 0.3] and z s S [0.9, 1.1]. 



prediction for that model. With a smaller statistical error 
and/or residual bias and cosmic variance, the true model 
with eg = 0.6 will become excluded from our 3 — a con- 
tours in the right panel of figure ^. In fact a maximum 
likelihood analysis on the noise free catalogue would give 
as = 0.8, that is 20% larger than the true as, which cor- 
responds to the lack of power in the predicted non-linear 
signal. We will get back to a more detailed discussion in 
Section 6 about this problem. 

5.2. VIRMOS-DESCART data 

We first consider flat cosmologies, since this is the class of 
models currently favore d by the cosmic microwave back- 
ground measurements ( de Bernardis et al., 2000 ), but 
alternative open universes are also investigated. In ei- 
ther case, the likelihood function is computed for 17640 
(10 x 14 x 14 x 9) models using Eq.(||), as a function of an- 
gular scale, and for a regular spacing in the default prior 
box. 

Figure ^ shows the four parameters constraints for dif- 
ferent priors and marginalized parameters for the flat cos- 
mology. The dashed lines shows the 68%, 95% and 99.9% 
contours when the default prior is applied for the two re- 
maining parameters. We cannot extract strong constraints 
in that case, but the right panel shows an interesting ten- 
dency between T and z s : a flat power spectrum (large 
r) can account for an underestimated source redshift. 
The thick solid curves shows the same contours with a 
stronger prior: T and z s are marginalized over [0.1, 0.4] and 
[0.39, 0.54] for the left panel, and VLm and as are marginal- 
ized over [0.1,0.4] and [0.7,1.3] for the right panel. We 
obtain the following constraint from the left panel: 



4xl0- 5 



'2x10- 




10 20 
0[arcminj 



30 



Fig. 7. The filled circles with error bars show the aper- 
ture mass measured on the rCDM simulated sky images, 
while the filled small squares show the signal measured 
in the input catalogue. The open circles show the mea- 
sured residual B-mode. Measurements and simulation are 
in perfect agreement, but the non-linear prediction ob- 
tained from Peacock fc Dodds (1996] ) for this model (solid 
line) is significantly off. 



for the 68% level and 

as = (0.58 ± 0.13) ^^0.025)^-0.48 



(16) 



&8 



(0.57 ±0.04) ^ 24 =F0.i8)n„- 



0.49 



(15) 



for the 95% contour. Constraints on the mass power spec- 
trum can be obtained from the right panel if one assume 
that photometric redshifts provide the exact redshift dis- 
tribution (which is given by z s = 0.44). In that case we 
have T G [0.12,0.38]. 
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Fig. 8. Left: constraints on f2 and erg for the flat cosmologies. The confidence levels are [68, 95, 99.9] from the brightest 
to the darkest area. The galaxy sample has a magnitude cut mi > 21. The gray area, and the dashed contours 
correspond to the contours computed with a full marginalization over the default prior T G [0.05, 0.7] and z s G 
[0.24,0.64]. The thick solid line contours are obtained from the prior T G [0.1,0.4] and z s G [0.39,0.54] (which is a 
mean redshift z s G [0.8, 1.1]). Right: constraints on T and 17m for the flat cosmologies. The contours have the same 
statistical meaning as for the left panel, but here, the dashed lines correspond to a marginalization over the default 
prior n G [0.1, 1] and a s G [0.3, 1.6], and the thick lines for n G [0.1,0.4] and cr 8 G [0.7, 1.3]. 




1 10 

0[arcmin] 

Fig. 9. The aperture mass statistic (Mf p ) measured on 
the data (see Figure |l|) compared to all the models in- 
cluded in the 68% contour (shaded area). For each mea- 
surement point, the error bars from left to right are: statis- 
tical errors, statistical error and residual bias, statistical 
errors and bias and cosmic variance. 

Figure ^| shows the aperture mass measurements with 
all the models inside the 68% contours as the shaded area. 
The error bars show the contribution the three errors as 



This work 68% 
Pierpaoli et al. 68% 
Seljok 68% 
Viana et al. 95% 
Pen 68% 




Fig. 10. The f2— as constraints for a flat universe from our 
work, compared to the cluster normalization constraints. 



variance added. We see that the statistical noise dominates 
at small scale, while the systematic residuals dominate at 
larger scales, the cosmic variance is never an important 
contribution. 

It is interesting at this stage to compare our results 
with measurements from other surveys. A comparison 
with Cosmic Microwave Background (CMB) constraints 
reveals that weak lensing will be helpfull to break the 
degen eracy between a& and Ojvf. Recently, Lahav et al 



a funct ion of scale. Each set of errors shows 3 bars, which (2001 ) have shown CMB estimation of these two param 



from left to right are: statistical noise, bias added, cosmic eters, assuming that the Hubble constant is a Gaussian 
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Fig. 12. Same as Figure H for the open cosmologies. 



0.10 0.20 0.30 0.40 0.50 0.60 0.70 

r 




Fig. 11. Solid line: the CMB alone constraints as de- 



scribed in Lahav et al. (2001 ) (figure borrowed from their 
paper) 



The CMB priors arc given by a concordance 
model, ' with the Hubble constant h — 0.7 with an r.m.s. of 
0.07, the baryon density fif, = 0.02, the primordial spec- 
tral in dex n — 1, and the rconization depth r — 0. Dashed 



lines o ur cosmic 3hcar constraints for the flat cosmologies 

(see Figure [|). In either case, the contours show the 68% 
and 95% confidence levels. 



variable centered at h = 0.7 with an r.m.s. of 0.07, and 
fixing other parameters (primordial spectral index n = 1, 
baryon density = 0.02 and reionization depth r = 0). 
Their results are shown as the solid lines on Figure [Tl]. 
An overlay of our constraints on the same plot (dashed 
lines) show that a combination of CMB and lensing would 
favorise low density models (Qm ~ 0.3 — 0.4) and rather 
low normalisation (as ~ 0.7 — 0.8). This plot reveals that 
the lensing constraints are a lmost orthogonal to the CMB 
constraints. As explained in Lahav et al. (2001 ), a weaker 
prior on h would extend the CMB contours and restore 
the degeneracy between Qm and erg, making Qm = 1 a 
viable solution again. But we see that lensing rules out 



such a solution because Qm ~ 1 with as > 0.8 is ex- 
cluded. Given that CMB alone predicts a flat Universe, 
the inconsistency between CMB and lensing for SIm = 1 
should be interpreted as in favor of a non-zero cosmologi- 
cal constant. The fact that CMB and lensing have oposite 
constraints in the (cts^m) parameter space make them 
indeed very complementary. A complete analysis which 
take into account the marginalization over the other pa- 
rameters (baryon density, r, etc..) is under way. However, 
we should note the agreement between our results and 



the combined CMB+2dF contraints (Lahav et al., 2001 
Figure 5). 

We should also compare our results more closely to the 
cluster normalization constraints, since those two methods 
are expected to probe a similar combination of as and £Im ■ 
Figure Uu shows our results and those obtained from clus 



tor measurements. As it was claimed before (Maoli et al 



2001; Van Waerbeke et al., 2001b; Rhodes et al., 2001), a 



joint estimate of Qm and as from weak lensing is consis- 



tent with the former cluster abundance estimates (Pen 



1998; Pierpaoli et al., 2001). Recently, these estimates 



were revisited but the new results are puzzling (see Figure 
0): whereas [Beljak (200l"l ) is marginaly consistent with 
our constraints, on the other hand, Viana et al. (200~i| ) 
is significantly lower. Pen (199"s| ) performed direct hydro- 
dynamic simulations to predict the cluster X-ray temper- 
ature function for various cosmological models. This by- 
passes the difficult mass ladder, of converting N-body or 
Prcss-Schechter mass functions into a temperature func- 
tion, and/or accounting for scatters in the relation, and 
the results are in good agreement with the cosmic shear 
constraints. However, some effects that may still not be ac- 
counted for in simulations include non-gravitational feed- 
back from galaxies, magnetic fields, thermal conduction, 
which may all limit the intrinsic accuracy of cluster nor- 
malizations. We will not enter into the debate between the 
cluster estimates here, but if the low normalization is con- 
firmed, this discrepancy might be an important finding: it 
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might be an indication of the inaccuracy of the non-linear 
predictions, as shown in Section 5.1. 

The maximum likelihood analysis was also carried out 
for open cosmologies. The probability contours shown in 
Figure ^| summarize the results which are indeed similar 
to the flat case. However, low density (SIm < 0.2), open 



interpretation is supported by ray-tracing simulations in a 
tCDM model, and more generaly by a comparison of the 
VIRGO simulations with the Peacock & Dodds non-linear 
prescription. It was claimed a 15% accuracy in their origi- 
nal paper ( Peacock fc Dodds, 1996j ), although it might be 



a bit more for some cosmological models (Jenkins et al 



universes, seem more difficult to reconcile with the data 1998). This is clearly the maximum uncertainty we can 



than flat models. This contradiction between observations 
and low density open universe results mostly from the 
small scale measurement of (M^ p ) (Figure^). Indeed, low 
density open universes predict too much power at small 
scale as compared to what can be allowed from the am- 
plitude of (M% p ) on scale of about one arc-minute. This 
clear difference between op en and flat ACDM un iverses 
was already pointed out by Schneider et al. (1998| ) (Fig. 
3) but this is the first time that it manifests on real data. 

Finally, although uncertainties are still large enough to 
leave room for a large sample of models, it is interesting to 
show how cosmic shear data can be used jointly with sev- 
eral independent surveys (not only CMB). Because weak 
lensing analyses probe dark matter in a direct way, cosmic 
shear are the best suited surveys to constrain as ■ It is then 
relevant to only focus on this parameter, using values of 
other cosmological parameters as they are derived from ex- 
ternal data sets. Assuming the mean rcdshift of sources is 
z s = 0.9 and h = 0.7±0.1 ( |Frccdmanct al., 200l| , from the 
HST Key Project), a flat universe (from CMB data) with 
a baryon fraction inferred from BBN and V = 0.2 ± 0.05 
( [Bzalay et al., 2002] , from the SDSS redshift survey), we 
then have Q M ~ 0.3 ± 0.1. In that case, the VIRMOS- 
DESCART cosmic shear survey provide a$ = 0.98 ± 0.06, 
in good agreement with other independent methods. As 
compared to the join CMB-cosmic shear alone discussed 
previously, the normalisation is higher. This is mainly due 
to the low Qm (i.e. high erg) combined with a strong prior 
on r (inferred from the galaxy redshift survey) and on the 
source redshift. 

6. Conclusion 

We explored a 4-dimensional parameter space using the 
most recent cosmic shear data. We included all possi- 
ble sources of error: statistical noise, cosmic variance and 
residual systematics. We obtained constraints on 17m 5 the 
power spectrum slope T, its normalization ag and the 
redshift of the sources z s . We marginalized over T and 
z s . Both the marginalization, and the inclusion of all 
the sources of error, make our results for (flM,cr&) ro- 
bust. We pointed out the complementarity between cosmic 
shear and CMB measurements for breaking the degener- 
acy among J7m and as , and the good agreement with CMB 
and CMB+2dF constraints. However, our results are only 
in marginal agreement with the latest cluster abundance 
constraints, which give a lower normalization a% ~ 0.7 for 
Ojvf ~ 0.3. If this discrepancy is confirmed in either mea- 
surements, this could be interpreted as an indication that 
the lensing non-linear prediction is not accurate enough 
given the already small size of cosmic shear errors. This 



tolerate with today's lensing measurements, and it will be 
insufficient for forthcoming surveys. This potential prob- 
lem suggests three paths for improvements: 

1. Ray-tracing simulations should be used more inten- 
sively to test non-linear schemes. So far, simple ray- 
tracing have been performed, by putting all the sources 
at a single redshift (z s = 1), although in principle we 
should expect a redshift dependence of the non-linear 
predictions failure (because different physical scales 
are probes for a varying redshift, for a fixed angular 
scale) . 

2. Progress on the theory side should be done. There 
might be some hope by reviving the halo models which 
give predictions close, but not identical, to Peacock 
& Dodds (like the peakpatch approach, Bond, priv. 
comm.). 

3. Ultimately, cosmic shear observations will lead to a 
measurement of the 3D mass power spectrum in non- 
parametric way, and therefore solve all the problems 
associated with non-linear modeling. This is possible 
only if the cosmological parameters are determined 
by other means. For instance, the linear mass power 
spectrum and cosmological parameters measurement 
at large scales using combined (or not) lensing data 
with cosmic microwave background, X-rays, could be 
obtained, and used to deconvolve the non-linear power 
spectrum. This means that we will be able to decon- 
volve the projected mass power spectrum measured 
from cosmic shear observations and recover the true 3D 
power spectrum. This is a work in progress, in which 
we are trying to recover the galaxy-galaxy and galaxy- 
mass corre lations as well, using tomography techniques 
( |Hu, 1999|) . 
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